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Executive  Summary 

The  five  manuscripts  reproduced  in  this  progress  report 
cover  the  research  of  the  Principal  Investigator,  Dr.  Dan 
Shechtman,  and  his  collaborators  under  the  DARPA/ONR 
sponsored  research  program  on  icosahedral  quazicrystals 
formed  by  rapid  solidification  of  alloys.  First, 

-microscopic  evidence  for  quasi-periodicity  in  an  alloy  with 
long-range  icosahedral  order  is  presented.  Local  5-fold 
rotational  axes  are  observed  to  be  uniformly  distributed 
down  to  the  atomic  scale.  Twinning  and  modulated  structures 
are  ruled  out  based  on  the  experimental  data. 

r In  a  high  rjasolution  electron  microscopy  study  on 

Mar.;  ,  f,'  o  pi  .t»G_ 

rapidly  sol idi f ied 'AlgMn  alloy  ribbons  the  relevant 
characteristics  of  quasi-periodicity  were  demonstrated,  The 
particular  topological  properties  of  such  a  periodic  network 
are  best  observed  by  direct  imaging  of  the  quasi-lattice  in 
the  electron  microscope. 

The  indexing  of  icosahedral  quasiperiodic  crystals  is 

x 

described  by  Cahn,  Shechtman  and  Gratias  including  the  step 
necessary  to  prove,  by  diffraction,  that  an  object  is 
quasiperiodic.  Various  coordinate  systems  are  discussed  and 
an  explanation  is  given  for  choosing  one  aligned  with  a  set 
of  three  orthogonal  two-fold  axes.  Using  this  coordinate 
syste:;  the  main  crystallographic  projections  are  presented 
and  several  analyzed  single-crystal  electron  Jiff  action 
patterns  are  demonstrated. 
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In  another  manuscript  Shechtman  provides  experimental 
evidence  for  quasiperiodic  crystals^  The  experiments 
include  various  diffraction  techniques  for  studying  the 
long-range  order  as  well  as  the  methods  to  determine  the 
local  atomic  order.  The  diffraction  pattern  has  well  defined 
sharp  peaks  and  fivefold  rotational  symmetry  and  the 
crystals,  in  differeht  orientations,  have  more  of  these  five 
fold  diffraction  patterns.  The  lattice  imaging  technique 
was  used  extensively  to  study  lattice  defects.  The 
technique  can  be  used  to  detect  the  fine  structure  of 
dislocation  cores  and  microtwin  boundaries.  No  boundaries 
were  seen  4n  the  samples  under  investigation  while  the 
quasipefr iodic  sequence  of  planes  is  clearly  visible. 

Neutron  diffraction  was  used  by  Kozer,  Cahn,  Gratias 
and  Shechtman  to  investigate  the  icosahedral  phase  of 
aluminum-manganese  alloys.  All  the  peaks  appear  at  the 
angles  in  agreement  with  the  icosahedral  indexing  with  a 
six-dimensional  cubic  lattice  parameter  of  approximately 
0.65  mm.  that  decreased  with  increasing  Mn  content. 

The  research  results  furnish  additional  strong 
documentation  for  the  existence  of  the  icosahedral 
quasicrystals.  A  variety  of  characterization  techniques  are 
described  and  the  experimental  data  is  analyzed  and 
explained.  The  work  advances  our  knowledge  and 
understanding  of  quasicrystals  and  their  relationship  to 
chemical  composition  and  solidification  conditions. 


MIlROSujFI  EVIDENCE  FOP  QUASI -F'ER  I  GDI  C  IT'i  IN  A  SOLID 
UITH  LONG-RANGE  I COSAHEDRAL  ORDER. 

Abstract : 

Higl'i  resolution  electron  microscopy  images  of  an 
icosahedral  solid  reueal  many  of  the  important  geometrical 
characteristics  of  quasi -periodici  ty  .  In  addition,  the 
local  ^i— fold  rotational  axes  are  seen  to  be  uniformly 
distributed  down  to  the  atomic  scale;  we  find  no  evidence  of 
an  underlying  3~ dimensi onal  lattice,  thereby  ruling  out  the 
hypotheses  of  microtwinning  and  modulated  structures. 
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introduction  of  six  indices,  two  per  independent  unit 

vectors  ,  the  fiert  one  associated  w  1  t  h  integer  length,  the 

f 

second  one  with  the  irrational  length  corresponding  to  the 
golde  mean.  Z  (£=1.6130  34...),  an  algebraic  irrational 

number  equal  to  2cosJL/5  and  therefore  directly  induced  by 
the  fiue-fold  symmetry  (10).  A  discrete  Fourier  spectrum 
with  a  finite  number  of  integer  indices  but  greater  than  the 
dimensionality  defines  a  structure  tc<  be  quasi -per  i  odi  c(  3 )  . 
An  important  property  of  the  quasi— periodicity  results  in 
the  unusual  progressing  of  the  intensities  of  the 

reflections.  More  a-' d  more  intense  reflections  occur  at 

distances  n+mT  in  which  m/n  are  the  successive  convergent 
approxirnants  of  L  .  The  geometrical  properties  of 

quasi  per  i  odi  c  patterns  have  extensively  studied  these  last 
years  (for  a  general  survey  of  aperiodic  tilings  see 
reference  11);  rigourous  algebraic  derivations  are  now 
proposed  (12,13,14;  which  generalize  the  original  Penrose 
tilings  (  4 )  ,  All  patterns  or  five  —  fold  s  y  rnm  etry  exhibit 
essential  characteristics  like  Fibonacci  sequences  and 
invariance  through  homothetic  multiplication  by  the  golden 
mean  . 

The  high  resolution  image  has  indeed  all  the  essential 
properties  of  a  two-dimensional  generalized  Penrose  tiling: 
short  Fibonacci  sequences  of  two  fundamental  lengths  related 
by  the  golden  mean  are  identified  both  in  the  spacinq  of 
rows  and  in  the  sp  aci n  g  of  bright  spots  within  rows;  also 


the  5  srfiB  sequences  and  pattern?  occur  at  eeoeral  different 
scale;  o»  -ngth.  This  is  particularily  cleat  for  pentaqone 
of  br 1 gh  spots  which  interpenetrate  at  larger  and  larger 
scales.  Although  there  is  no  strict  periodicity,  any  finite 
feature  on  the  micrograph  is  found  repeated  and  oriented  the 
same  way  at  some  finite  (but  not  constant!)  distance.  This 
property-  constitutes  the  essence  of  qua si- periodic  patterns 
which  are  the"  most  periodic"  among  the  aperiodic  patterns. 
We  therefore  conclude  that  we  arc  seeing  in  Figure  1  a 
projection  along  a  five-fold  axis  of  a  quasi -per i odi c 
3-dimensional  structure  with  icosahedral  symmetry.  This 
conclusion  is  consistent  with  the  recent  illuminating 
mathematical  derivation  of  Katz  and  Duneau  (14)  who 
demonstrated  that  the  cuts  along  all  high  symmetry 
orientations  (two-,  three-  and  five-fold  axes)  of  an 
icosahedral  quasi-lattice  are  indeed  2-dimensional 
quasi-lattices.  A  complete  report  of  the  irf-aging  of  all 
high  symmetry  zones  in  high  resolution  microscopy  will  be 
published  elsewhere  (15). 

The  structural  homogeneity  of  the  icosahedral  phase- 
down  to  the  atomic  scale  definitely  contradicts  the 
hypothesis  of  micrctwi  nmng  in  the  sense  in  which  this 
concept  is  commonly  used.  If  rrucrotuinninq  is  redefined  to 
include  twinning  of  a  size  below  the  resolution  of  this 
microscope,  the  concept  will  have  acquired  a  new  meaning, 
which  would  raise  more  questions  than  it  would  solve.  It 
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would  still  hs'ie  to  e-' plain  the  long  range  orientational 
order,  the  typical  Fibonacci  sequences,  the  spatial 

f 

self-simi lar 1 ty ,  as  well  as  the  unusual  diffraction 
intensities  which  are  perfectly  well  reproduced  in 
kinernatical  calculations  based  on  quasi  -per  1  odi  ci  ty 
(13,14,15)  . 

Meak  beam  electron  microscopy  does  reveal  a  mottling  in 
the  5nm  range.  This  has  been  interpreted  as  evidence  for 
the  presence  of  mi  crocr  yst  als  and  hence  of  twinning  (8); 
from  our  high  resolution  study  ue  attribute  this  contrast 
mottling  to  the  local  variations  of  strains;  this  is 
confirmed  by  a  continuous  tilting  experiment  in  the 
microscope:  the  bright  area  in  dark  field  images  move 
smoothly  with  no  discontinuity. 

A  more  subtle  problem  is  the  distinction  between  this 
structure  and  i ncommen sura tel y  modulated  structures.  Ne 
found  no  evidence  of  an  underlying  lattice  in  direct  space 
which  is  consistent,  in  the  reciprocal  space,  with  the  fact 
that  the  diffraction  pattern  can  not  be  partitioned  into 
"main"  and  "satellite"  reflections.  As  far  as  the 
mathematical  aspects  are  concerned,  any  quas i -per i odi  c 
function,  including  one  with  icosahedral  symmetry,  may  be 
regarded  as  the  diagonal  cut  of  a  purely  periodic  function 
of  N  variables  where  N  is  the  dimension  of  the  algebraic 
base  (here  N  =  3x2)  .  Accordingly,  Janner  and  Janssen  (16) 


the  case  for  orientations  of  the  N-d  i  men  si  o  rial  space  which 
correspond  to  isolated  strata:  as  a  result  the 
i  ncommensur  at  i  on  is  locked  to  a  fixed  cut  angle  and  cannot 
vary  with  temperature  and  composition.  This  is  the  case  for 
the  icosahedral  phase  which,  in  fact,  corresponds  to  a 
special  cut  of  a  6-dimensi onal  simple  cubic  (10,14).  These 
special  properties  set  the  icosahedral  phase  into  a 
particular  class  of  modulated  structures,  the  one  that 
corresporids  to  isolated  strata  whose  little  groups  do  not 
belong  to  the  32  crystallographic  point  groups:  the 
icosahedral  phase  is  a  "special  point"  long  range  ordered 
phase  with  noncryst allogr aphi c  little  group.  This  analysis 
suggests  that  the  icosahedral  phase  could  be  a  STABLE  phase 
under  certain  conditions  of  composition  temperature  and 
pressure  exactly  like  usual  crystalline  phases. 
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Figure  1:  (a)  Electron  diffraction  of  the  5-fold 

(1 T  0)  zone  axis. 

(b)  Cor r espondi ng  high  resolution  image. 

(c) Optical  Fourier  transform  of  (b)  showing 
the  restoration  of  the  five-fold  symmetry. 
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High  Resolution  Electron  Microscopy 


of  Aluminum-based  Icosahedral  Quasi -Crystals 
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Introduction . 

Some  rapidly  solidified  binary  and  ternary  alloys 
exhibit  a  long  range  ordered  structure  with  no  translational 
periodicity  (Shechtman  et  al  1984,  1985, Zhang  et  al  to 
appear)  but  with  a  discrete  Fourier  spectrum  characteristic 
of  Almost-periodici ty  (Besicovith  1932).  These  have  been 
recently  called  "quasi-crystals"  (  Levine  et  al  1985)  and 
can  be  described  by  the  Cut  and  Projection  Method  (C.P.M. ) 
(Duneau  et  al  1985,  Elser  to  appear,  Kalugin  et  al  1985) 
which  is  a  generalized  discrete  version  of  the  earlier 
hyperspace  description  used  for  continuous  density  functions 
of  incommensurate  structures  (Janssen  et  al  1984).  The 
particular  topological  properties  of  such  aperiodic  networks 
are  best  observed  by  direct  imaging  of  the  quasi-lattice  in 
the  electron  microscope.  It  is  shown  here  that  the  relevant 
characteristics  of  quasi-periodicity  are  present  in  the 
images  which  can  be  interpreted  independently  from  the 
actual  organization  of  the  atomic  species  within  the 
quasi-per iodic  framework. 


Experimental . 

The  observations  reported  here  were  made  on  rapidly 
solidified  ribbons  of  the  A16Mn  alloy.  Electron  microscopy 
was  performed  on  a  JEOL  200CX  provided  with  the  Cs=lmm  pole 
piece  and  the  -+10°  goniometer  stage.  The  HR  images  along 
a  five-fold  axis  are  shown  on  Fig.l  for  different  defocus 
values  and  objective  apertures  .Although  the  images  show  no 
rigourous  five-fold  symmetry  (only  the  Patterson  of  the 
structure  does  show  an  exact  icosahedral  symmetry)  numerous 
homogeneously  distributed  pentagons  and  decagons  are 
observed  at  atomic  scale.  White  dots  in  Fig.l-a  fit 
temarkably  well  with  the  simple  projection  of  the 
quasi-lattice  generated  by  C.P.M. ;  this  is  confirmed  by 
Fig. 2  which  corresponds  to  the  very  important  two-fold 
orientation . 
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Figure  1:  (a)  5-fold  orientation  at  -90  Nm  defocus 
with  0.08  Nm  objective  aperture.  (b)  -170Nm  defocus,  same 
aperture . 
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Figure  1:  (c)  5-fold  orientation  at  -150  Nm  defocus 

with  0.03  Nm  objective  aperture. 

Figure  2:  2-fold  orientation  at  -50  Nm  defocus  with 
0.05  Nm  aperture. 
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Here  too, the  projected  quasi- latt ice  nodes  and 
intensity  modulations  of  the  images  correspond:  neither 
dynamical  effects  nor  phase-changes  due  to  aberrations  can 
modify  the  intrinsic  topology  of  the  images.  The  unit 
length  of  the  edge  of  the  elementary  rhorabohedra  (  Mackay 
1982,  Kramer  et  al  1984)  constituting  the  basic  tiles  for 
aperiodic  space  filling  is  .46  nm;  this  parameter  gives  a 
unit  length  of  .65  nm  for  the  6-dimensionnal  primitive 
hypercubic  generating  lattice.  These  values  strongly 
suggest  that  the  actual  structure  contains  certainly  more 
than  one  atom  per  unit  cell;  not  only  some  of  the  Al-Mn 
distances  in  the  usual  crystalline  phases  are  known  to  be 
abnormally  short  (Cooper  et  al  1966)  but  the  intrinsic 
volumes  of  the  rhombohedra  are  far  too  large  to  be 

consistent  with  the  experimental  density  (Kelton  et  al  1985) 
if  occupied  by  a  single  atom  (an  average  of  4  atoms  is 
expected).  From  the  strict  mathematical  point  of  view,  the 
C.P.M.  is  not  able  to  make  any  distinction  between  any  two 
homothetic  quasi-lattices  which  are  in  the  ratio  X.  5  (where 
T>  is  the  golden  mean).  This  property  is  directly  observed 
in  both  diffractions  and  images  of  the  two-fold  orientation 
/see  Portier  et  al  1985  for  a  detailed  discussion/  whereas 
an  apparent  x,  scaling  is  observed  on  the  five-fold 

orientation . 

Convergent  Beam  Patterns  of  the  different  principal 
orientations  all  show  an  uniform  contrast  within  the  discs. 
This  effect  might  be  explained  by  the  fact  that  the 

reciprocal  quas i- latt ice  being  a  dense  Z-modulus  (see 
below),  each  intensity  within  the  discs  results  from  an 

infinite  number  of  contributions  of  diffracted  beams.  An 
additional  plausible  explanation  is  the  fact  that  actual 
quasi -crystals  have  a  short  correlation  length  (Bancel  et  al 
1985)  due  to  an  imperfect  quasi-periodicity  in  the  material 
which  would  weaken  the  effect  of  the  excitation  parameters 
in  the  dynamical  diffraction. 


Discussion . 


The  very  challenging  problem  which  has  to  be  solved  is 
the  description  of  the  atomic  structure.  Electron 
microscopy  is  considered  as  being  a  suitable  tool  for 
collecting  informations  about  the  positions  of  the  atoms  and 
some  models  (Hiraga  et  al  1985,  Guyot  et  al  1985,  Knowles  et 
al  1985)  have  already  been  proposed  based  mostly  on  only  the 
five-fold  orientation.  There  however,  serious 
difficulties  in  interpreting  images:  the  dynamical 
calculations  depend  crucially  e  number  of  allowed 
diffracted  beams  (Cornier  et  ai  issue)  which,  contrary 
to  the  case  of  crystals,  form  a  dense  set.  The  so-called 
"small  divisor"  problem  (Belissard  1982)  in  the  almost 
periodic  diffraction  Hamiltonian  makes  the  perturbation 
expansion  diverge,  leading  to  approximate  solutions  which 
are  critically  dependent  on  the  chosen  cut-off  (Cornier  et 
al  this  issue)  .  For  an  optimal  cut-off,  the  simulated 
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dynamical  images  of  a  simple  quasi-lattice  with  an  average 
atom  at  each  quas i -lat t ice  node  fit  remarkably  well  the 
experimental  ones  (Cornier  et  al  this  issue).  Hence,  HR 
images  of  relatively  close  packed  structures  are  essentially 
governed  by  the  topology  of  the  lattice  and  are  quite 
insensitive  to  the  atomic  motif.  Indeed,  a  specific  result 
of  the  C.P.M.  is  that  the  Fourier  components  of  the 
potential  are  multiplied  by  the  Fourier  transform  of  a 
generic  cut  function  which  does  not  depend  on  the  atomic 
species.  This  term  is  so  predominant  that  chemically 
different  icosahedral  phases  -like  those  in  Shechtman  et  al 
1985  ,  Zhang  et  al  1985-  do  show  very  similar  intensity 
distributions  although  they  are  constituted  of  different 
atoms.  These  reasons  show  that  HR  microscopy  will  be  of 
limited  use  for  the  structure  determination  of  metallic 
quasi-crystals . 
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Since  the  definition  of  quasiperiodicity  is  intimately  connected  to  the  indexing  of  a  Fourier 
transform,  for  the  case  of  an  icosahedral  solid,  the  step  necessary  to  prove,  using  diffraction,  that 
an  object  is  quasiperiodic,  is  described.  V arious  coordinate  systems  are  discussed  and  reasons  are 
given  for  choosing  one  aligned  with  a  set  of  three  orthogonal  two-fold  axes.  Based  on  this 
coordinate  system,  the  main  crystallographic  projections  are  presented  and  several  analyzed 
single-crystal  electron  diffraction  patterns  are  demonstrated.  The  extinction  rules  for  three  of 
the  five  icosahedral  Bravais  quasilattices  are  compared,  and  some  simple  relationships  with  the 
six-dimensional  cut  and  projection  crystallography  are  derived.  This  analysis  leads  to  a  simple 
application  for  indexing  powder  diffraction  patterns. 


I.  INTRODUCTION  II.  QUASIPERIODICITY 

The  recent  discovery  of  Shechtmanite, 1,2  a  metallic  The  Fourier  transform  of  a  periodic  function  is  a  set 

phase  with  long-range  icosahedral  orientational  sym-  of  delta  functions  that  are  periodically  spaced  and,  in 

metry  and  experimentally  discrete  diffraction  patterns,  general,  vary  in  magnitude.  Diffraction  gives  informa- 

has  revealed  a  new  class  of  ordered  structures.  The  ico-  tion  in  the  form  of  a  Fourier  transform  of  the  correla- 

sahedral  symmetry  is  inconsistent  with  strict  crystallo-  tions  of  an  object.  If  the  object  is  a  periodic  crystal,  the 

graphic  periodicity,  yet  discrete  diffraction  implies  qua-  diffraction  pattern  is  a  discrete  set  of  spots  of  varying 

siperiodicity.  Schechtmanite  is  thus  cited  as  an  example  intensity  that  are  positioned  on  a  reciprocal  lattice, 

of  a  quasiperiodic  crystal  or  quasicrystal,  for  short.  Three  reciprocal  lattice  vectors  form  a  basis  to  locate 

Since  the  two  icosahedral  groups  are  not  part  of  the  32  any  spot  in  a  three-dimensional  reciprocal  space, 

crystallographic  point  groups,  the  possibility  exists  that  A  mathematical  function  is  quasiperiodic  by  defini- 

any  of  the  infinity  of  noncrystallographic  point  groups  tion  if  its  Fourier  transform  is  a  set  of  delta  functions 

will  be  observed.  Indeed,  claims  that  specimens  exhibit-  that  are  not  uniformly  spaced  as  they  would  be  for  a 

ing  two  other  point  groups,  decagonal  U/m  (Ref.  3)  periodic  function,  but  whose  spacing  can  be  described 

and  duodecagonal  12  (Ref.  4),  have  since  been  report-  by  a  finite  set  of  lengths. 5  If  an  infinite  number  of  lengths 

ed.  are  required  the  function  is  cdled  almost  periodic. 

These  quasicrystalline  phases  present  challenging  Therefore  specimens  that  give  countable  diffraction 

problems  in  crystallography.  In  this  paper  the  mundane  spots  that  cannot  be  indexed  with  three  reciprocal  lat- 

housekeeping  problems  of  coordinate  systems,  index-  tice  vectors  are  quasiperiodic  if  they  can  be  indexed  with 

ing,  and  extinction  rules  that  are  the  essential  language  a  finite  set. 

of  reporting  experimental  observations  will  be  dis-  It  has  been  shown  that  any  ^-dimensional  quasiper- 

cussed.  Indexing  is  not  just  a  housekeeping  procedure:  it  iodic  function  requiring  a  basis  of  N  vectors  can  be  con- 

is  an  essential  part  of  proving  that  a  structure  is  periodic,  sidered  to  be  derived  from  a  periodic  Ar-dimensional 

quasiperiodic,  or  almost  periodic.  We  will  concentrate  function  cut  by  a  D-dimensional  plane.3  If  every  spot  in 

on  the  icosahedral  phases  and  the  three-dimensional  the  diffraction  pattern  can  be  indexed  using  a  combina- 

aspects  of  descriptions  that  are  most  easily  derived  in  tion  of  N  reciprocal  lattice  vectors,  then  the  object  that 

higher  dimensions.  gave  this  diffraction  pattern  can  be  represented  by  a  cut 

of  a  iV-dimensional  periodic  object. 

‘ 1  Present  address:  Center  for  Materials  Research,  The  Johns  Hopkins  Thc  -COSahedral  point  group  IS  not  consistent  With 

University,  Baltimore,  MD  21218;  also  a  guest  worker  at  the  Na-  translational  periodicity.  The  icosahedral  diffraction 

tiona!  B'treau  of  Standards,  Gaithersburg,  MD  20899.  pattern  cannot  consist  of  periodically  spaced  spots.  If 
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the  diffraction  pattern  consists  of  discrete  spots  that  can 
be  indexed  by  a  finite  number  of  basis  vectors  it  is  quasi- 
periodic. 

Several  papers  have  already  discussed  ways  of  de¬ 
riving  quasipenodic  structures  from  cuts  of  periodic 
higher  dimensional  structures6-10  including  icosahedral 
structures.  The  converse  problem  of  taking  a  particular 
diffraction  pattern  and  indexmg  it  has  been  attempted 1 1 
in  a  way  that  has  been  criticized.12,13  The  prof  .m 
centers  on  the  fact  that  with  a  combination  of  inc  ora- 
mensurate  lengths  any  spot  can  be  located  approximate¬ 
ly  with  any  desired  degree  of  accuracy.  We  will  show 
that  with  our  indexing,  the  observed  high-intensity 
spots  form  a  simple  sequence  in  which  none  are  missing 
and  none  left  out. 


III.  COORDINATE  SYSTEMS 

Taking  basis  vectors  along  important  symmetry  di¬ 
rections  simplifies  the  crystallographic  formulation. 
For  any  group  with  a  unique  rotation  axis,  the  z  axis  is 
taken  parallel  to  that  axis.  This  is  the  proper  choice  for 
tne  crystallographic  groups  such  as  hexagonal  and  the 
noncrystallographic  groups  such  as  decagonal.  The  ico¬ 
sahedral  groups  have  6  fivefold  axes,  10  threefold  axes, 
and  15  twofold  axes.  Taking  one  of  the  fivefold  axes  as 
the  z  axis  (Fig.  1 )  leaves  the  other  five  in  a  ring  63.43° 
from  this  axis.  Although  these  all  would  make  acute 
angles  with  the  z  axis,  there  are  obtuse  angles  between 
some  of  them.  There  is  no  choice  of  sign  for  the  six  axes 
that  would  give  equal  angles  between  all  of  them.  A  co- 


0/1  1/0  1/1 


0/1  1/0  1/1 

FIG.  1  The  stereozraphic  projection  of  the  principal  symmetry  directions  and  mirror  planes  of  the  icosahedral  group  m35  as  seen  along  the 
fivefold  direction.  The  number  indicates  the  indexing  system  that  will  be  described  later  in  the  text. 
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ordinate  system  based  on  these  axes  is  not  only  skewed  The  simplest  system  is  a  cubic  coordinate  system  in 

but  encounters  the  difficulty  of  keeping  track  of  obtuse  which  the  axes  of  the  coordinate  system  are  aligned  with 

and  acute  angles.  There  are  right  angles  between  a  five-  a  set  of  three  orthogonal  twofold  axes  of  the  icosahedral 

fold  axis  and  five  of  the  twofold  axes,  that  could  serve  as  group  (Fig.  3).  The  15  twofold  axes  fall  in  five  such 

a  coordinate  system  when  comparing  the  icosahedral  all  equivalent  to  each  other  through  the  operation  of  a 

phase  with  the  decagonal  phase.  fivefold  axis.  This  is  the  coordinate  system  used  in  the 

If  we  examine  the  icosahedral  group  with  a  three-  International  Tables  of  Crystallography,  and  it  is  the 

fold  axis  along  z  ( Fig.  2 ) ,  we  have  three  choices  of  coor-  coordinate  system  we  will  use,  even  tc  .  escribe  the  other 

dinate  systems  that  are  related  to  familiar  crystallogra-  coordinate  systems.  It  has  all  advantages  of  orthogonal 

phic  ones.  The  6  fivefold  axes  now  fall  into  two  groups,  axes. 

either  of  which  could  be  used  as  a  rhombohedral  basis.  The  three  coordinate  systems  discussed  here  are  all 

In  one  set  the  fivefold  axes  make  acute  angles  with  each  based  on  subgroups  of  the  icosahedral  group,  as  shown 

other.  In  the  other  set  the  angles  are  obtuse.  A  hexagon-  in  Fig.  4.  Using  a  coordinate  system  based  on  a  lower 

al  coordinate  system  could  be  based  on  the  threefold  symmetry  than  icosahedral,  requires  special  attention 

axis  and  the  three  twofold  axes  at  right  angles  to  it.  Of  for  the  icosahedral  symmetries  not  used  in  that  coordi- 

special  interest  are  the  three  twofold  axes  at  71°  from  z  nate  system:  equivalent  reflections  will  not  necessarily 

that  form  an  orthogonal  set.  have  similar  indices. 


0/1  1/0  1/1 


3-fold  0/1  1/0  ~u~ 


FIG,  2.  The  stereographic  projection  of  the  icosahedral  group  m 35  as  seen  along  a  threefold  direction.  Note  the  possibility  of  a  rhombohedral 
coordinate  system  using  either  set  of  fivefold  axes,  a  hexagonal  coordinate  system  using  the  twofold  axes  along  the  equator  or  the  three  twofold 
are*  at  71*  The  thickened  great  circles  show  the  three  zone  axes  of  the  diffraction  pattern  indexed  in  Figs.  5-7,  namely  the  twofold  labeled 
(0/0  1/00/1 ),  the  fivefold  labeled  [0/0  1/0  0/1  j  and  the  threefold  axis  labeled  [0/0  1/2  0/1]. 
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IV.  THE  CUBIC  COORDINATE  SYSTEM 

Having  selected  a  coordinate  system  aligned  along 
three  perpendicular  twofold  axes,  we  need  to  describe 
vectors  and  planes  in  both  the  direct  and  reciprocal 
spaces.  An  immediate  advantage  of  the  choice  of  a  cubic 
coordinate  system  is  that  the  indices  of  a  piane  and  the 
direction  normal  to  it  are  the  same.  We  define  unit 
lengths  a  and  a*  in  both  spaces  and  dimensionless 
lengths  and  position  in  terms  of  these.  Consider  a  posi¬ 
tion  ( UVW)  or  a  plane  (HKL).  The  set  of  all  positions 
or  planes  equivalent  through  the  operations  of  the  icosa¬ 
hedral  groups  (235)  and  m35  is  as  follows: 

1.  The  threefold  axes  give  cyclic  permutation,  but 
note  that  while  ( UVW),  (  WUV),  and  (  VWU)  are  equi¬ 
valent,  (  VUW),  etc.,  is  different. 

2.  The  twofold  axes  give  pairs  of  sign  changes,  the 
mirrors  give_  individual  sign  changes.  Thus  for  235, 
(UVW),  (UVW),  ( UVW),  and  (UVW)  are  equivalent 
but  ( UVW),  etc.,  is  not.  For  m35  all  sign  changes  are 
equivalent  positions. 

3.  The  fivefold  axes  introduce  a  change  of  the  mag¬ 
nitude  of  U,  V,  and  W and  introduce  the  golden  section 
t  —  2  cos  36°  =  ( 1  +  v 5 ) / 2  =  1.618  034.  Because  the 
International  Tables  have  a  misprint  in  the  rotation  ma¬ 
trices,  we  repeat  them  here  in  Table  I,  using  their  nota¬ 
tion  G  =  r/2,  g  =  1/2 r  =  G  —  1/2.  Performing  the  ma¬ 
trix  multiplication  we  obtain 

(HKL)Y 


=  \\(H-K)+r(K-L)t 
(L  —  H)  +  t(H  +  K), 

-  (K  +  L)  +  r(L  +  H) ).  (1) 

In  this  multiplication  we  make  frequent  use  of  the  iden¬ 
tities 


It  may  at  first  seem  surprising  that  equivalent  posi¬ 
tions  turn  out  to  have  designations  with  different  nu 
merical  values  of  the  components,  but  this  is  unavoid 
able  with  a  Cartesian  coordinate  system  and  groups 
with  such  a  high  symmetry.  The  most  symmetric  crys¬ 
tallographic  group  mlm  has  order  48.  It  has  48  unit 
triangles.  Taking  all  the  permutations  and  sign  changes 
of  (UVW)  generates  48  equivalent  general  positions, 
one  in  each  unit  triangle.  It  is  impossible  to  generate  the 
60  or  120  equivalent  positions  needed  (resp.  for  235  or 
m 35)  with  just  three  symbols.  Thus,  the  choice  of  a  cu¬ 
bic  coordinate  system  results  in  the  possibility  that  as 
many  as  five  different  sets  of  indices  may  be  necessary  to 
represent  equivalent  positions  or  planes. 

Now  consider  a  plane  or  reciprocal  space  position 
with  an  index  of  the  form 

(h  +  h  'r,k  +  k  'r,l  +  I't), 

that  is,  where 

H-h-Yh'r,  K  =  k  +  k V,  L=l+  I't, 

in  which  the  h,  h',  k,  k  ' ,  l,  and  / '  are  all  integers.  We 
introduce  the  six-index  notation  (h/h',  k/k',  l //') 
or  (h  /h  '  k  /k'  l // ')  to  designate  such  a  reflection. 

Operation  of  the  fivefold  rotation  will  change  the 
numerical  value  of  the  six  integers. 

(h  +  h  'r,k  +  k  'r,l  +  l'r)Y 
=  \((h  +  k')  -  (*  +  /') 

+  r[(  -/+*')  +  (k  —  /')], 

(  -h  +k‘)  +  (/  +  h’) 

+  T\(h  +  k  )  +  (£  +  /’)], 

(-*  +  /')  +  (/  +  *') 

+  r[(h-k')  +  (/+/>')]). 


(5) 


t2  =  1  +  r, 

1/r  =  r  —  1, 
G2  -l-g2  +  (j)2 


(2) 

(3) 

(4) 


In  order  for  this  to  be  a  close  set,  we  impose  the 
parity  conditions  that  the  three  sums  h  +  k',k  + l',  and 
/  +  h  'are  even:  after  fivefold  rotation  the  six  integer  in¬ 
dices  remain  integers  and  the  parity  rules  are  conserved. 


TABLE  I.  Matrices  for  fivefold  rotation  about  (1  r  0] . 


fi 

g 

C\ 

/ 

*  G 

1  \ 

Y=Y  <=(* 

G 

-l 

r3  =  k-3  = 

G  1 

-g  ) 

\G 

J 

*/ 

\ 

-  j  * 

-o 

~g 

E5=y-2=  G 


-g 


r  = 


G  =  (v5  +  1  )/4  =  r/2  =  cos  36’  =  0.809017 
g=(fi-  D/4  =  \/2t  =  G-  j  =cos  72”  =  0.309017 
Kr=  Y  ' 
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1  k 

Tf 


Finding  that  all  spots  in  the  diffraction  pattern  are  of 
this  form  implies  quasiperiodicity. 

We  next  turn  to  the  indexing  of  the  principal  direc¬ 
tions  ( Figs.  1-3 ) .  The  six  fivefold  axes  are  all  of  the  form 
( IrO),  or  in  the  six-index  notation  (1/00/1  0/0). 

The  ten  threefold  axes  have  two  different  designa¬ 
tions:  four  are  along  (111)  or  (1/0  1/0  1/0)  and  six 
along  (r2  1  0)  or  ( 1/1  1/00/0).  Note  that  both  vectors 
have  been  chosen  to  have  a  length  ^3. 

Three  of  the  twofold  axes  are  along  the  cube  axes. 
Theremaining  12  have  the  form  ((7,g,I)  or  (0/1  T/l  1/ 
0). 

V.  RECIPROCAL  QUASILATTICES 

Consider  a  quasilattice  in  three-dimensional  reci¬ 
procal  space  in  which  every  spot  occurs  as  a  sum  of 
integer  multiples  of  a  finite  number  ( greater  than  three ) 
of  vectors.  We  will  compare  two  lattices  formed  only  of 
equivalent  vectors.  In  particular,  let  us  first  take  the  six 
vectors  along  the  fivefold  axes,14 

Q  =  2 /ml,  (6) 

i  =-  l 

where  the  n,  are  integers  and  the  q,  are  ( IrO)  which  in 
the  six-index  notation  is 
Ql  =  (1/0  0/1  0/0), 
q2  =  (0/1  0/0  1/0), 
q3=  (0/0  1/0  0/1), 
q4  =  (T/0  0/1  0/0), 

q,=  (0/10/0  1/0), 

q6  =  (0/0  I/O  0/1). 

[There  are  384  ways  for  choosing  the  q , ’s.  All  give  equi¬ 
valent  results.  The  choice  of  the  set  (7)  corresponds  to 
th  two  rhombohedral  bases,  q,,  q2,  q3  define  the  acute 
rhombohedron  and  q4,  q,,  q6  the  obtuse  one  (see  Fig. 
2).  ]  The  set  of  six  numbers  ( n , )  can  be  considered  an 
indexing  of  Q,  and  has  been  used  in  a  number  of  papers. 
This  same  set  can  also  be  considered  to  be  a  six-dimen¬ 
sional  lattice  vector.  To  express  Q  in  terms  of  the  three- 
dimensional  cubic  coordinates  we  substitute  set  (7)  in 
Eq.  (6)  and  perform  the  summation, 

Q  =  ((«,  -n4)/(n2+nf),  (n3  -n6)/(ni  +  n4), 

(n2  —  ns)/(n3  +  n6)).  (8) 


We  can  thus  convert  from  the  six -dimensional  vec¬ 
tor  ( n ■  n2,  n3,  n4,  nf,  n6)  to  the  six-index  three-dimen 
sional  vector  (h/h',k/k',l /l'). 

h  =  rii  —  n4,  h'  =  n2  +  n}, 

k  —  n2  —  n6,  A: '  =  n,  +  n4,  (9a) 

l  =  n2-n5,  1'  —  n2  +  n6, 
and  vice  versa 

2nl  =  h  +  k',  2n4=  —h  +  k', 

2n2  =  l  +  h',  2n5  =  —  /  +  k  ',  (9b) 

2n3  =  *  +  /',  2  n6=  -k  +  r. 

The  form  of  Eq.  (9b)  demonstrates  that  Q  in  the  six- 
index  notation  obeys  the  parity  rules.  Ihese  restrictions 
on  the  indices  are  extinction  rules  and  are  given  in  Table 
II  under  the  heading  of  the  P  quasilattice  ( to  be  defined 
later). 

These  parity  mles  lead  to  four  kinds  of  positions: 

1.  six  even  indices, 

2.  four  even  indices  (odd/even,  even/odd,  even/ 
even), 

3.  two  even  indices  ( even/odd,  odd/even,  odd/ 
odd), 

4.  six  odd  indices. 

If  the  coordinates  of  a  spot  (HKL )  that  is  consistent 
with  the  parity  rules  are  multiplied  by  r,  the  results  are 

(tHtKtL)  =(h'/(h  +  h'),k  /(k  +  k'), 

/'/(/  +  /')).  (10) 

Only  if  the  original  (HKL)  =  {h  /h'  k/k'  l /l') 
contains  two  or  six  even  indices,  will  (  tH  tK  tL  )  sa¬ 
tisfy  the  parity  rules.  On  the  other  hand,  scaling  by  r3 

(r'Hr'Kr3!.)  =  ((A  +  2h  ')/(2h  +  3/t '), 

(k  +  2k')/{2k  +  lk'), 

(/  +  2/')/(2/+3/'))  (11) 

preserves  the  parity  rules  for  all  {HKL)  that  satisfy 

them.  These  scaling  rules  are  tabulated  for  the  primitive 
(P)  reciprocal  quasilattice  in  Table  III. 

In  addition,  the  square  of  Q  is  of  the  form 

02  =  A'  +  ;Vfr,  (12) 

where  using  Eqs.  ( 2 )  and  ( 3 ) , 
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TABLE  III.  Parity  and  scaling  properties  for  the  P  and  F  reciprocal 
quasila  trices. 


Parity  of 
indexes 

Scaling 

P 

F 

Q1  =  N+tM 

N  M 

six  even 

r 

Present 

Present 

4  n 

4  m 

four  even 

■r3 

Present 

Absent 

4  n  +  2 

4m  +  1 

two  even 

r 

Present 

Present 

4  n 

4  m 

six  odd 

r3 

Present 

Absent 

4a  +  2 

4m  +  1 

N  =  lj 

'  n)  =  h  2 

* 

+  h 12 

+  k2  +  k 

'2  +  l2  +  I 

’2,  (13) 

i  -  i 


M  =  h'2  +  k 12  + 1'2  +  2(hh'  +  kk'  +11').  (14) 

Equation  (13)  indicates  that  N,  the  integer  part  of  Q 2,  is 
twice  the  square  of  the  length  of  the  six-dimensional  n 
vector  and  therefore  must  be  even.  Because  of  the  parity 
rules,  if  (Vis  divisible  by  4,  M  also  is  divisible  by  4,  and  if 
N  is  not  divisible  by  4,  M  is  of  the  form  4m  -f  1.  In 
addition,  as  we  will  show  later,  M  conforms  to  the  limits 
—  N/t  <M<Nt  (15) 

and  that  the  most  intense  reflections  occur  for  the  lar¬ 
gest  value  of  M,  which  we  shall  denote  by  M0.  For  N 
divisible  by  4 

M0  =  4  L  (Nt)/4  J  ,  (16) 

where  L  x  J  is  the  largest  integer  in  x 
When  N  is  not  divisible  by  4, 

M  =  4m  -f  1  <  Nt. 

The  largest  value  of  the  integer  m  is 
m0  =  L  (Nt  —  l)/4  J 

and 

M0=  1  +4L((Vr—  l)/4  1 .  (17) 

We  next  define  a  two-parameter  indexing  Q(N,r) 


Q2(N,r)=N  +  M0T-4rr, 

r  —  0,1,...  <  L  (N/t  +  M0)M  J  ,  (18) 

in  which  the  Q0(N)  will  turn  out  to  define  the  sequence 
of  intense  reflections. 

QI(N)=Q2(N,0)=N  +  Mo(N)t.  (19) 

The  same  procedure  has  been  used  by  us  to  find  the 
quasilattice  formed  from  the  15  vectors  along  the  two¬ 
fold  axes.  The  same  result  is  obtained  in  a  simr  ’«r  man¬ 
ner  by  putting  the  face-centering  restrictions  f  on  the 
six-dimensionai  lattice  formed  by  the  n, ,  and  then  using 
Eq.  (9a) .  The  result  for  that  and  for  the  body-centered  I 
lattice  are  also  given  in  Table  II. 

In  Table  III  we  make  a  simple  comparison  between 
the  P  and  F  which  share  the  same  parity  rules.  It  can  be 
seen  that  the  F  lattice  scales  uy  r,  and  has  only  N  —  4n, 
the  P  lattice  comprises,  in  addition,  N  =  4n  +  2  spots 
that  scale  by  r3.  Figure  5  is  a  diffraction  intensity  calcu¬ 
lation  using  the  sphere  approximation  in  the  cut  and 
projection  method.8,9 

VI.  INDEXING  THE  SINGLE  QUASICRYSTAL 
PATTERNS 

Figures  6-8  show  the  indexing  of  the  electron  dif¬ 
fraction  patterns  of  Refs.  1  and  2.  The  stereographic 
projection  (Fig.  2)  presents  the  three  zone  axes  of  these 
diffraction  patterns,  namely,  the  [0/0  0/0  0/2]  two¬ 
fold,  the  [0/0  1/0  0/1  ]  fivefold,  and  the  [0/0  1/2  0/1  ] 
threefold  axes.  Note  that  all  of  the  reflections  in  the  five¬ 
fold  and  threefold  zone  axis  conform  to  either  the  For  F 
reciprocal  quasilattices  of  Sec.  V.  These  could  not  be 
used  to  distinguish  between  the  two  reciprocal  quasilat¬ 
tices. 

It  is  easy  to  show  that  the  reflections  for  which 
h  '  +  k  '  +  /'  is  odd  v  i>!  appear  in  neither  the  five-  or 
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FIG.  5.  Diffraction  patterns  calculated 
by  the  cut  and  projection  method  for  all 
three  reciprocal  cubic  quasilattices  with 
the  same  lattice  parameter  a*  for  the 
two-  and  fivefold  zone  axes. 
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FIG.  8.  Indexed  diffraction  pattern  of  the  [  0/0  0/0  0/2 1  zone  axis.  Note  the  square  array  made  of  all  even  and  two-even  rpots  that  scale  with  r.  If 
only  these  spots  were  present,  this  pattern  would  have  fourfold  symmetry  and  scale  by  r. 
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threefold  zone  axes.  For  example,  in  order  for  a  reflec¬ 
tion  to  be  in  the  [0/0  1/00/1]  zone  axis  K  4-  tL  =  0. 
This  implies  that  £  +  /'  +  r(  &'  +  /  +  /')  =  0,  which 
can  only  occur  if  (&  +  /')  =  (k'  + 1  +  /')  =  0.  Adding 
the  even  number  h'  —  l  gives  h  '  +  k '  +  l '  —  In.  The 
same  result  is  obtained  with  the  [  1 1 1  ]  zone  axis.  Fitting 
the  two  zone  axes  alone  can  not  distinguish  between  P 
and  F.  Indeed  the  model  of  Levine  and  Steinhardt  seems 
to  fit  only  the  three-  or  fivefold.14  In  the  [0/0  0/0  0/2] 
zone  axes  /  =  / '  =  0,  but  h  +  k  ( and  h'  4-  k ' )  can  be 
odd.  This  gives  the  spots  with  two  odd  indices,  that  can 
be  seen  in  Fig.  7.  Spots  with  six  odd  indices  derive  from 
these  after  some  fivefold  rotations.  The  twofold  zone 
axes  thus  show  some  spots  with  r3  scaling. 

The  twofold  zone  axis  for  the  F  quasilattice  show 
only  spots  with  r  scaling  and  may  show  an  accidental 
fourfold  symmetry.  Permutation  of  the  x  and  y  indices 
preserves  the  parity  rules  and  does  not  change  Qe  ( see 
Sec.  VII).  If  the  intensity  is  only  a  function  of  Qc  the 
pattern  will  show  the  fourfold  axis. 

Observing  an  icosahedral  diffraction  pattern  in 
which  all  the  spots  can  be  indexed  with  this  six  index 
notation  proves  that  we  have  a  diffraction  pattern  of  a 
quasipenodic  structure.  In  particular  the  experimental 
pattern  is  not  the  For  I  Furthermore,  because  it  can  be 


indexed  with  six  integers,  the  object  can  be  represented 
as  a  slice  and  projection  of  a  six -dimensional  periodic 
structure. 

Several  other  methods  of  indexing  with  six  numbers 
have  been  successful."-13  They  have  differed  from  each 
other  in  ways  that  become  significant  after  we  examine 
the  cci  and  projection  method. 


VII.  THE  CUT  AND  PROJECTION  FROM  SIX 
DIMENSIONS 

In  three  dimensions  the  six  vectors  along  the  five¬ 
fold  axes  are  not  orthogonal.  We  can  choose  the  six¬ 
dimensional  cubic  space  in  which  each  of  these  vectors  is 
a  basis  vector  along  a  hypercube  axis  that  is  perpendicu¬ 
lar  to  all  the  others.  The  set  of  the  six  numbers  nt  then 
represents  a  position  vector  in  the  six-dimensional  cubic 
reciprocal  lattice,  and  the  equations  (8)  give  the  corre¬ 
spondence  between  positions  in  the  three-  and  six-di- 
mensional  spaces.  The  cut  and  projection  is  accom¬ 
plished  by  rotating  the  six-dimensional  space  so  that 
what  will  become  three  axes  in  the  three-dimensional 
space  are  in  the  cut  plane.  The  six-dimensional  rotation 
matrix  corresponding  to  the  chosen  set  (7)  is 
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1  T  0-1  T  0“ 

T  0  1  T  0—1 

0  1  T  0-1  T 

T  1  0  T  1  0 

1  0  -T  1  Or 

0  —  T  1  0  1  r_ 

(20) 

Tlie  three-by-six  matrix  constituting  the  top  half  of 
the  rotation  matrix  consists  of  the  six  fivefold  directions 
in  three  dimensions  as  column  vectors  [compare  with 
Eq.  ( 7 )  ] ,  or  the  coordinates  of  the  three-dimensional 
cube  axes  in  the  six-dimensional  space  as  row  vectors. 
The  bottom  three  by  six  matrix  consists  of  the  comple¬ 
mentary  orthogonal  space  lost  by  the  cut.  It  consists  of 
projection  directions  from  six  dimensions  into  the  real 
space.  In  the  six-dimensional  reciprocal  space  the  dis¬ 
tance  of  a  spot  from  the  cut  plane  is  related  to  its  intensi¬ 
ty  in  the  three-dimensional  reciprocal  space. 

There  are  five  Bravais  lattices  consistent  with  icosa- 
hedral  symmetry  in  six  dimensions. 10  Three  of  these  are 
the  P,  I,  F  hypercubic.  Consider  first  the  primitive  lat¬ 
tice.  Its  reciprocal  iattice  is  also  primitive.  Because  the 
projection  of  each  six-dimensional  lattice  vector  is  a 
(1/00/1  0/0)  vector  along  the  fivefold  axes,  the  reci¬ 
procal  quasilattice  observed  experimentally  to  be  one 
composed  of  such  vectors  will  be  denoted  P.  The  other 
lattices  in  our  example  correspond  to  the  six-dimension¬ 
al  face-centered  F  and  body-centered  /  reciprocal  lat¬ 
tices  F corresponding  resp.  to  /  and  F  direct  lattices. 

In  the  P  six-dimensional  lattice  the  location  of  a 
spot  in  the  dimensionless  units  that  we  have  used  in 
three  dimensions  is,  using  Eq.  ( 20 ) : 

Q6  =  yj  2(1  +  r2)(rt,,rt2,«3.«4,rt5,rt6).  (21 ) 


FIG  9.  Two-dimensional  representation  of  the  cut  and  projection 
procedure  for  the  magnitude  of  Q  and  the  distance  from  the  cut  plane 

Q 


As  a  result,  using  Eqs.  ( 2 )  and  (13), 

Ql  =  2(2  +r)(/if  +n\  +  n\  +  n\  +  n]  +  n\) 

=  (2  +  t)N.  (22) 

The  length  of  Q6  is  related  to  the  N  part  of  the  length  of  Q 
(Fig.  9).  Since  projection  shortens  Q 

Q2<Ql-  (23) 

Using  Eqs.  (12)  and  (22) 

N  +  Mr  <  (2  +  t)N. 

Therefore 

M<  [(1  -f  t)/t]N  =  tN.  (24) 

Since  Q 2  >  0,  we  also  have 

—  N /t<M  (25) 

proving  the  inequality  ( 15).  If  we  define  Qc  to  be  the 
distance  of  a  spot  in  the  six-dimensional  reciprocal  space 
from  the  cut  plane  (Fig.  8),  we  have 

Q]  =  Ql -Q2 

=  t(Nt  —  M).  (26) 

Thus  the  largest  M  for  a  given  N  will  have  the  smallest 
Qc- 

Up  to  now  we  have  been  using  a  dimensionless  Q. 
We  introduce  a  three-dimensional  quasilattice  constant 
d0  such  that  the  three-dimensional  diffraction  vector  k 
and  the  interplanar  spacing  are  given  by 

k  =  Q/d0,  (27) 

d(h  /h  ',  k/k  ',1/1')  =d0/^N+rM .  (28) 

The  length  of  the  six-dimensional  reciprocal  lattice  con¬ 
stant  a?  is  related  to  d0  by 

a*d0  =  \  2(2  +  t)  .  (29) 

In  Table  IV  we  give  the  values  of  Q0(N)  and  the 
corresponding  Qc  that  according  to  theory  should  be 
inversely  correlated  with  intensity.  This  sequence  of  Q ’s 
has  an  obvious  beginning  (N  —  2)  and  using 
d0  —  1.7466  nm  we  produce  for  0. 155  nm  radiation  two 
further  columns,  the  diffraction  vector  k  and  the  diffrac¬ 
tion  angle  20.  Comparing  these  with  our  own  and  pub¬ 
lished  powder  diffraction  data11  indicates  a  one-to-one 
correspondence  with  this  list.  The  only  omissions  are 
due  to  overlap  with  fee  aluminum  and  low  intensities  in 
the  published  data. 1 1  The  justification  for  our  choice  of 
d0  is  to  match  Table  IV  with  experimental  data.  The 
choice  of  d0  in  Ref.  1 1  was  different  and  this  leads  to 
complications  to  be  discussed. 

In  Table  V  we  list  the  Q0(N)  series  in  three  different 
notations  each  using  six  indices.  The  («,■)  and 
{h  /h'  k  /k'  l  // ' )  are  representative  of  that  particular 
Qq(N).  The  six-dimension  (n, )  vector  is  that  which 
projects  into  the  longest  possible  three-dimensional  vec- 
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TABLE  IV.  Scattering  variables  for  the  P  lattice. 


Dimensionless 

Dimensional 

(</„  =  1.74657) 

N 

V0  Ql 

Qc 

d  '  ( nm  ' ) 
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TABLE  V.  Vanous  indexing  methods  for  the  first  12  strong  reflections. 


2 

1 

(100000) 

4 

4 

(100100; 

6 

9 

(111000! 

8 

12 

( 101 101 ) 

10 

13 

(llioil; 

12 

16 

(210010) 

(mill) 

14 

21 

(201101) 

16 

24 

(211011) 

18 

29 

(211111) 

20 

32 

(201201) 

22 

33 

(211201) 

24 

36 

(222000) 

(211211) 

(1/0  0/1  0/0) 

12 

(2lll 1 1 ) 

(0/0  0/2  0/0) 

30 

(220011) 

(1/1  1/1  1/1) 

20 

(110001) 

(321112) 

(0/0  2/2  0/0) 

30 

(UlOlO) 

(1/2  2/10/0) 

60 

(221020) 

(2/2  0/2  0/0) 

60 

(311111) 

0/2  2/2  0/0) 

12 

( 1/0  2/3  0/0) 

60 

(21T001) 

(331021) 

(2/2  2/2  0/0) 

60 

(211101) 

( 1/2  2/3  0/0) 

12 

(100000) 

[ ( 32T002 )? ] 

(0/0  2/4  0/0) 

30 

(110000) 

(0/1  2/4  1/0) 

120 

(2/2  2/2  2/2) 

20 

(220002) 

(0/2  2/4  0/0) 

60 

[(561033)?] 

tor,  and  the  two  are  related  by  Eq.  (9).  For  powder 
pattern  intensities  the  multiplicity  of  each  spot  is  given. 
The  sequences  of  our  indexing  of  course  satisfies  Eq. 
(13)  relating  jV  to  the  sum  of  squares  of  indices.  As  was 
noted  earlier,  the  list1 1  of  powder  diffraction  angles  cor¬ 
relates  perfectly  with  this  sequence,  with  only  one  omis¬ 
sion  in  this  range. 

The  last  column  gives  the  indexing  of  Ref.  1 1 .  The 
fundamental  ( 100000)  vector  was  chosen  in  Ref.  1 1  to 
have  the  length  of  Q0(  1 8 ) .  As  a  result,  all  smaller  reflec¬ 
tions  become  higher  index  and  the  simple  monotonic 
relation  between  N  and  Q  of  Eq.  ( 1 3 )  is  lost.  Because  of 
this  there  is  no  obvious  start  to  the  sequence,  and  it  is 
difficult  to  know  if  any  intense  spots  are  missing.  Fur¬ 
thermore,  because  of  the  incommensurability  one  can 
approach  any  angle  with  arbitrary  precision  by  using 
high  indices.  The  two  assignments11  labelled  with  ques¬ 
tion  marks  are  probably  such  approximations.  The 
question  of  which  vector  one  chooses  as  a  fundamental 
length  in  this  case  is  uniquely  resolved  by  the  sequence 
of  the  intense  reflections.  Choosing  a  longer  vector  leads 
to  the  problems  cited  above.  It  will  be  difficult  to  choose 
a  shorter  vector,  because  intensities  are  likely  to  be  very 
low. 

We  have  so  far  concentrated  on  the  Q0  series.  Table 
V  lists  representative  spots  for  all  the  spots  that  occur 
for  N  =  12  spherical  shell  in  six  dimensions  (Fig.  9). 
Two  types  of  indices  sum  to  N  =  12.  In  six  dimensions 
there  are  26  spots  of  type  (111111)  and  236!/(2!3!)  of 
type  (211000).  Six  different  lengths  result  from  project¬ 
ing  in  three  dimensions  ranging  from  Ql  =  37. 89  to  Q , 
=  5.53.  The  longer  Q's  are  from  six-dimensional  vec¬ 
tors  nearly  parallel  to  the  cut  plane;  the  shorter  ones  are 
from  those  nearly  perpendicular.  By  themselves  the 
( n, )  give  no  clue  about  projected  length  until  the  projec¬ 


tion  has  been  defined  in  Eq.  (20),  but  in  the  cubic  in¬ 
dexes  the  lengths  obviously  become  shorter  as  the  2’s 
shift  from  primed  to  unprtmed  positions  and  sign  differ¬ 
ences  appear  in  the  indexes. 

The  ( 1 1 )  1 1 1 )  are  high  symmetry  axes  in  six  dimen¬ 
sions.  They  project  into  three  dimensions  as  either 
three-  or  fivefold  axes  with  four  different  lengths.  The 
(211000)  project  onto  mirror  planes  or  (110)  planes 
with  six  different  lengths.  In  powders  the  Q0  ( 12)  reflec¬ 
tion  will  be  a  superposition  of  72  individual  diffraction 
spots.  These  multiplicities  are  most  readily  apparent  by 
forming  ratios  of  the  cubic  indices  and  comparing  these 
with  the  indexes  of  the  three  symmetry  axes  and.  the 
mirror  planes. 

We  hope  to  have  demonstrated  that  the  cubic  coor¬ 
dinate  system  has  many  advantages  over  either  the  skew 
coordinate  system  or  the  six-dimensional  one.  In  addi¬ 
tion,  there  is  an  obvious  simplification  when  a  particular 
lattice  constant  is  chosen. 

VIII.  DISCUSSION 

Several  indexing  methods  have  been  introduced 
and  need  to  be  compared.  We  have  introduced  a  method 
of  indexing  based  on  a  three-dimensional  cubic  coordi¬ 
nate  system  using  icosahedral  symmetry.  Six  indexes  are 
necessary'  and  sufficient,  indicating  that  the  icosahedral 
solid  is  quasiperiodic  and  can  be  represented  as  an  irra¬ 
tional  cut  of  a  six-dimensional  periodic  solid.  Much  of 
the  geometry  is  developed  without  recourse  to  six  di¬ 
mensions. 

The  six-index  notation  is  merely  a  shorthand  for 
indexing  irrational  numbers  of  the  form  h  +  h  t.  Ordi¬ 
nary  vector  addition  generates  the  three-dimensional 
diffraction  pattern  from  a  single  vector  replicated  by  the 
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operations  of  the  icosahedral  group.  There  is  a  single 
basic  length:  the  ir/5  rotations  introduce  the  algebraic 
number  r  and,  therefore,  incommensurability.  Unlike 
other  incommensurate  structures,  in  which  the  incom¬ 
mensurate  ratio  varies  with  temperature  and  composi¬ 
tion,  the  ratio  r  is  part  of  fivefold  rotation  and  remains 
the  same  for  large  changes  in  composition15  and  lattice 
parameter  and  a  wide  variety  of  alloy  systems.  The  cubic 
coordinate  system  has  many  advantages  that  derive 
from  an  orthogonal  coordinate  system.  It  incorporates 
the  cubic  subgroup  part  of  the  icosahedral  symmetry, 
but  it  leaves  the  problem  of  multiple  notation  for  equiva¬ 
lent  positions. 

The  multiple  notation  is  really  not  a  problem.  We 
have  found  two  convenient  methods.  For  a  particular 
index,  multiplication  by  the  rotation  matrices  quickly 
gives  all  the  other  notations  for  equivalent  reflections. 
When  dealing  with  a  large  number  of  indexed  positions, 
they  are  quickly  sorted  by  calculating  Q 2  and  sorting  by 
iVandAf  (or  r).  Just  as  (330)  and  (41 1 )  coincide  in  bcc 
powder  patterns,  this  sorting  will  put  nonequivalent 
spots  in  the  same  (N,r)  box.  However,  since  all  the  spots 
in  the  same  ( N,r )  box  originate  from  the  spots  in  the  six¬ 
dimensional  reciprocal  space  that  are  the  same  distance 
from  the  cut  plane,  they  are  likely  to  have  similar  inten¬ 
sities.  The  other  three-dimensional  coordinate  systems 
have  skewed  axes  aligned  along  equivalent  symmetry 
axes.  They  have  all  the  problems  of  skewed  axes  includ¬ 
ing  the  problem  that  permutation  of  indices  not  only 
cause  length  changes  but  can  cause  the  vector  to  rotate 
to  a  nonequivalent  direction.  In  the  cubic  indexing,  be¬ 
cause  the  indices  are  primed  and  unprimed,  cyclic  per¬ 
mutations  of  ( HKL )  can  only  lead  to  equivalent  vec¬ 
tors.  Odd  permutations  do  not  change  the  length. 

In  six  dimensions  there  is  an  orthogonal  system,  but 
the  rotation  matrix  that  defines  the  cut  plane  requires 
arbitrary  choices  of  the  signs  of  the  fivefold  basis  vectors 
in  three  dimensions.  Thus  the  various  six-dimensional 
vectors  of  type  (110000)  represented  two  different 
lengths  in  three-dimensions  differing  by  a  factor  of  r, 
depending  on  the  relative  orientation  with  respect  to  the 

TABLE  VI.  A  complete  listing  of  the  .V  =  12  reflections. 


cut  plane,  or  equivalently  whether  the  angles  between 
the  various  ( 100000)  type  (010000)  vectors  in  three 
dimension  are  acute  or  obtuse. 

There  is  a  geometric  way  of  understanding  the  dis¬ 
crepancy  between  the  various  indexing  methods.  It 
might  seem  that  any  reflection  along  the  fivefold  axis 
could  serve  to  define  the  unit  length,  but  it  must  satisfy 
several  criteria.  The  first  is  that  all  other  spots  must  then 
be  indexable  with  six  integers  that  obey  the  parity  rules. 
The  spots  that  survive  this  criterion  differ  from  each 
other  in  length  by  powers  of  r 3.  In  six  dimensions  all  of 
these  vectors  lie  on  the  same  (xyyyyy)  fivefold  plane. 
The  shortest  distance  in  this  plane  is  ( 100000)  and  this 
is  the  one  which  must  be  found.  Only  one  of  the  fivefold 
reflections  corresponds  to  this  minimum  distance  and 
generates  an  orthogonal  basis  in  six  dimensions.  The  one 
chosen  in  Ref.  1 1  happens  to  be  the  (21 1 1  ll),  which 
leads  to  a  skewed  six-dimensional  basis.  This  destroys 
the  hypercubic  geometry  that  is  so  important  to  the  sim¬ 
plicity  of  the  indexing  we  have  proposed.  As  a  result  the 
indexing  fails  to  fit  the  criterion  of  a  simple  hierarchy. 
There  are  inifinitely  many  reflections  along  the  fivefold 
axis,  both  longer  and  shorter  than  the  one  which  meets 
this  criterion.  For  the  P  quasilattice,  it  is  the  longest 
fivefold  reflection  that  is  shorter  than  any  intense  reflec¬ 
tion  in  any  direction.  All  shorter  reflections  are  then 
projections  of  longer  six-dimensional  vectors  with  long 
Qc  and  weak  intensities. 

These  criteria  not  only  point  to  a  natural  and  unique 
indexing  but  solves  the  important  problem  of  what  con¬ 
stitutes  the  unit  reciprocal  lattice  vector.  Comparison  of 
Table  V  indicates  that  while  a  different  choice  of  this 
unit  still  gives  a  completely  consistent  set  of  indexing, 
there  is  a  simplicity  and  completeness  to  the  choice 
based  on  the  strong  reflections.  The  Q0  series  fits  the 
observed  intense  reflections  without  omission  and,  qual¬ 
itatively,  is  inversely  correlated  with  Qc ,  the  distance 
from  the  cut  plane.  The  next  in  the  sequence  (Table  VI) 
£?'(tV,  1)  =  Qo  (A0  —  all  have  Qc  >2 yr  which,  of 
course,  is  longer  than  Qc  for  any  of  the  Q0  spots,  and  are 
readily  distinguished  from  the  main  sequence  by  their 
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weak  intensity.  Our  indexing  method  therefore  uniquely 
identifies  the  quasilattice  parameter  d0  in  a  way  that  is 
quite  equivalent  to  that  chosen  by  Elser1*  based  on  an 
examination  of  the  Fourier  transform  of  the  cut  func¬ 
tion  along  a  systematic  row.  Because  our  method  is 
purely  geometric  and  does  not  refer  to  any  specific  cut 
function,  it  is  in  a  sense  more  general. 

In  this  paper  we  have  taken  for  granted  that  the 
symmetry  is  truly  icosahedral.  The  choice  of  the  cubic 
axes  for  a  coordinate  system  does  not  imply  that  we 
believe  that  this  is  a  cubic  crystal,  as  has  been  suggested 
by  others. 16,17  We  will  address  this  issue  in  a  separate 
paper.18 
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Abstract-  A  review  of  the  experimental  evidence  which  rejects 
the  determination  of  the  Icoaahedral  phase  as  periodic  is 
presented.  The  experiments  discussed  include  various  diffrac¬ 
tion  techniques  for  the  study  of  the  long  range  order  as  well 
the  methods  to  determine  the  local  atomic  arrangement. 


I  -  INTRODUCTION 

The  discovery  of  the  Icosahedral  phase  (1),  (2)  has  generated  a  great 

deal  of  research  activity.  It  is  the  fascination  of  being  a  part  of 
an  emerging  new  field.  which  motivated  investigators  from  several 
disciplines  to  contribute  about  150  articles  in  the  first  year.  How 
ever,  a  certain  amount  of  mainly  passive  scepticism  still  exists 
within  the  ranks  of  crystal lographers . 

Since  the  first  X-ray  diffraction  Experiment  by  von  Laue  in  1912,  all 
crystals  studied  by  various  techniques  could  be  defined  as  periodic. 
Many  structures  were  analysed  to  have  rather  complicated  motifs  and 
large  unit  cells,  which  contain  in  several  cases  well  over  a  thousand 
atoms,  (Frank-Kasper  phases  are  a  typical  example). 

All  these,  however,  could  be  shown  to  have  periodic  translational 
symmetry.  Therefore,  a  crystal  which  generates  sharp  diffraction 
peaks  was  axiomed  over  the  years  to  be  periodic.  This  axiom  which  has 
become  a  cornerstone  in  crystallography  is  not  supported  by  the 
mathematics  of  the  nature  of  diffraction.  On  the  contrary,  it  was 
demonstrated  and  proved  time  and  again  that  the  Fourier  transform  of 
almost  periodic  and  quasiperiodic  functions  also  generate  sharp 
peaks . 

The  question  whether  quasiperiodic  crystals  exist,  narrows  down  to  the 
question  whether  quasiperiodic  arrangement  of  atoms  can  form 
whether  it  can  be  stable. 

The  mathematical  as  well  as  other  theoretical  and  experimental  tools, 
to  deal  with  quasiperiodic  crystals,  were  available  for  quite  some 
time  before  the  announcement  of  the  discovery  of  the  Icosahedral 
Phase,  and  these  were  put  to  use  almost  immediately  following  the 
announcement . 
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The  experiments  that  revealed  the  exceptional  nature  of  the 
Icosahedra 1  phase  were  repeated  and  this  not  only  confirmed  the  first 
results,  but  also  expanded  them  and  au'ded  evidence  for  quas iper lodi- 
city  in  the  Icosahedral  phase  and  other  phases  (3)  (4).  The  phase, 
first  found  in  several  aluminum  alloys  (5)  was  later  discovered  in 
other  compositions  (6). 

The  purpose  of  this  article  is  to  discuss  experimental  evidence  which 
rejects  the  possibility  of  explaining  the  structure  of  the  Icosahedral 
phase  on  principles  of  classical  crystallography,  thus  creating  the 
foundation  for  the  understanding  of  the  structure  of  quasiper iodic 
crystals . 

EXPERIMENTAL 

There  is  no  doubt  that  electron_microscopy ,  with  its  versatility  and 
wealth  of  experimental  procedures,  played  a  key  role  in  the  early 
stages  of  the  study  of  the  Icosahedral  phase  and  other  quasiperiodic 
structures. 

The  property  that  drew  attention  to  the  fivefold  orientation  in  April 
1982  was  the  very  dark  image  of  the  Icosahedral  grains  which  were 
oriented  with  the  fivefold  axis  parallel  to  the  electron  beam.  The 
Icosahedral  phase  diffracts  electrons  in  an  unusual  way. 

Next  came  the  diffraction  pattern  (Fig.l).  with  well  defined  sharp 
peaks  and  fivefold  rotational  symmetry.  The  pattern  has  more  unique 
features,  the  distances  from  the  centrally  transmitted  beam  to  the 
diffracted  beams  are  related  by  various  powers  of  the  golden  mean  and 
the  intensity  of  the  diffracted  beams  does  not  decay  as  a  function  of 
the  distance  from  the  center.  But,  perhaps,  the  most  intricate 
observation  was  that  the  crystal,  in  different  orientations,  has  more 
such  fivefold  diffraction  patterns.  The  analysis  of  the  symmetries  of 
the  crystal  lead  to  the  composition  of  patterns  shown  in  Fig. 2  and  to 
the  realization  that  the  crystals  possess  Icosahedral  symmetry.  The 
explanation  of  the  previously  unobserved  set  of  symmetries  could  be 
either  by  using  the  knowledge  of  traditional  periodic  crystallography 
or  by  proposing  something  else  and  new. 


Fig.l.  Selected  area  diffraction  Pattern  taken  from  a  single 
Icosahedral  Crystal  (1). 
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Tig. 2.  The  sequence  of  selected  area  diffraction  pattersn  obtained 
from  one  Icoaahedral  crystal  at  various  tilt  angles  (1). 


Periodic  crystallography  may  suggest  that  each  spot  in  the  unique 
patterns  is  formed  by  a  contribution  of  several  differently  oriented 
particles  of  a  periodic  phase,  such  as  multiple  twinning.  To  examine 
this  idea  several  experiments  were  carried  out  in  the  electron 
microscope.  The  first  is  a  set  of  dark  field  images  obtained  from 
different  diffracted  beams  in  a  two-beam  condition  (1).  In  each  dark 
field  image  the  whole  Icoaahedral  grain  was  illuminated,  and  down  to 
the  resolution  of  a  few  nanometers,  no  twins  can  be  observed.  This 
experiment  was  repeated  since  on  various  other  Icoaahedral  phases, 
some  of  which  (such  as  several  rapidly  solidified  Al-Mn-Si  alloys!  are 
with  much  less  strain  contrast  with  the  same  results  (7). 

Convergent  beam  diffraction  patterns  were  taken  at  the  various 
important  orientations.  This  technique  involves  a  conical  electron 
beam  that  illuminates  a  spot  on  the  specimen,  the  diameter  of  which  is 
between  1.5nm  to  20  nm.  The  diffraction  patterns  obtained  this  way 
are  identical  to  the  ones  observed  by  the  selected  area  method,  and 
are  similar  to  one  another  across  the  Icoaahedral  crystal.  Such 
selected  area  diffraction  patterns  were  taken  at  various  thicknesses 
of  the  specimen,  and  down  to  a  few  nanometers  in  thickness,  indicating 
that  if,  indeed,  the  diffraction  pattern  is  a  composite  of  several 
patterns  obtained  from  several  periodic  crystals,  then  these  crystals 
must  be  very  small. 

Perhaps,  the  most  convincing  set  of  experiments  in  the  stuuy  of  the 
structure  of  the  Icoaahedral  phase  is  based  upon  lattice  imaging.  The 
lattice  imaging  technique  has  been  used  extensively  in  the  last  decade 
to  study  lattice  defects.  It  has  been  demonstrated  that  the  technique 
can  detect  the  fine  structure  of  dislocation  cores  and  that  of  many 
kinds  of  boundaries,  including  microtwin  boundaries  (8).  The 
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technique  was  applied  to  the  Icosahedral  phase  in  several  alloys  <9> 
revealing  the  atomic  structure  at  different  orientations.  Even  though 
it  is  difficult  to  confirm  a  structural  model  based  on  these 
observations,  several  facta  are  clear.  The  most  important  observation 
is  that  no  boundaries  can  be  seen  in  the  structure.  The  quasipenodic 
sequence  of  planes  is  clearly  seen,  and  the  orientations  of  planes  in 
the  2-,  3-  and  5-fold  orientation  is  conspicuous.  An  optical 
diffraction  taken  from  the  TEM  plate  to  form  a  Fourier  transform  of 
the  image  recreates  the  electron  diffraction  pattern.  This  result  is 
obtained  even  with  a  rather  small  aperture,  i.e.  with  the  diameter  of 
the  area  which  contributes  to  the  optical  diffraction  in  the  order  of 
several  nanometers.  The  various  rotational  symmetries,  including  the 
fivefold,  are  thus  a  fundamental  property  of  the  atomic  order  and  do 
not  result  from  any  composition  of  small  periodic  crystals,  such  as 
microtwins. 

X-ray  diffraction  has  traditionally  been  a  most  powerful  and  precise 
tool  in  the  study  of  crystals.  The  technique  was  put  to  use  in  the 
study  of  the  quasipenodic  crystals,  starting  with  the  Icosahedral 
phase  <1)(5)<10>.  The  fit  of  the  X-ray  diffraction  peaks  to  those  of 
the  electron  diffraction  was  confirmed,  and  the  indexing  of  the 
patterns  followed  soon  (5H11).  Even  though  there  is  no  dispute  as  to 
the  meaning  of  each  peak  among  the  various  suggestions  for  indexing 
languages,  the  indexing  techniques  vary  from  one  to  the  other.  The 
most  fundamental  of  all  is  the  one  suggested  by  Cahn  et  al  (11).  It 
is  based  upon  a  set  of  three  orthogonal  2-fold  axes  and  presents  an 
easy  and  workable  indexing  Bystem.  Another  important  result  of  the 
analysis  of  the  X-ray  diffraction  patterns  is  the  introduction  of  the 
step  necessary  to  prove,  using  diffraction,  that  a  solid  is  quasi¬ 
period  ic  . 

Among  the  other  results  obtained  by  X-ray  diffraction  is  the  accurate 
measurement  of  planar  spacing,  which  is  not  obtainable  accuiately  by 
electron  diffraction,  and  the  introduction  of  indexed  stereographic 
projections  for  the  Icosahedral  phase  (11).  In  general,  the  peaks 
observed  in  all  the  X-ray  studies  are  rather  wide  representing  a 
correlation  length  of  10-30nm.  Numerous  speculations  as  to  the  cause 
of  this  width  have  been  presented  including  local  strain,  faults  and 
uneven  distribution  of  atomic  constituents  to  name  a  few.  This 
question  is  unsettled  at  present  and  an  agreed  upon  explanation  is 
still  sought. 

Neutron  diffraction  on  powdered  Icosahedral  phase  was  also  performed 
on  three  alloy  compositions  (12).  Peak  positions  were  found  at 
identical  positions  to  those  previously  found  by  X-ray  diffraction  and 
previously  as  produced  by  cuts  and  projections  of  a  6D  cube  with  a 
lattice  parameter  of  about  0.65nm  (see  later  in  these  proceedings). 
The  intensities  of  the  neutron  diffraction  pattern  provide  information 
on  the  chemical  order  of  the  Icosahedral  phase  by  comparison  to  those 
of  the  X-ray  diffraction  pattern.  This  is  based  upon  the  difference 
between  the  structure  factors  characteristic  of  the  two  techniques, 
which  is  vividly  demonstrated  for  the  Al-Mn  system. 

The  local  atomic  structure  of  the  Icosahedral  phase  was  studied  by 
NCR,  XAFS  and  NMR. 

Shortly  following  the  discovery  of  the  Icosahedral  phase  Nuclear 
Gamma-Ray  Resonance  1.NGR2  experiments  were  performed  on  Al-  <Mn,Fe) 
alloys  (13).  The  NGR  spectrum  of  an  iron  atom  in  the  Icosahedral 
phase  is  a  probe  of  the  local  environment  of  that  atom  and  the 
technique  is  most  sensitive  to  the  first  one  or  two  near-neighbor 
shells.  A  detailed  analysis  of  the  spectra  obtained  from  the 
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Icosahed ra 1  structure  ghowa  that  they  can  be  fitted  to  a  pair  of 
symmetric  doublets,  indicating  two  non-symmetr ic  positions  for  the 
iron  atom.  It  is  assumed  that  the  iron  and  manganese  atoms  occupy  the 
small  site  types  and  are  distributed  arbitrarily  within  the  quasi- 
penodic  crystal.  The  results  do  not  suggest  a  model,  but  reject  any 
model  in  which  the  manganese  atom  occupies  a  symmetric  position. 

In  another  study  (14)  the  local  atomic  environment  in  the  Icosahedral 

phase  was  studied  by  Extended_X-ray_Absorpt ion  Fine_St ructure _ (EXAFS_>. 

measurements.  The  results  which  provide  information  that  lead  to  a 
structural  model,  demonstrate  the  power  of  this  tool  for  the  study  of 
the  order  of  the  Icosahedral  phase.  The  model,  detailed  in  these 
proceedings,  supports  the  randomly  connected  Icosahedra  model 
previously  suggested  (1)  and  suggests  a  structural  unit  which  consists 
of  a  cage  of  Mn  atoms  positioned  on  the  vertices  of  an  Icosahedran. 
The  model  also  suggests  that  the  quasiperiodic  crystals  grow  along  the 
threefold  axis,  which  is  consistent  with  previously  reported  optical 
microscopy  observations  (10). 

Nuclear  Magnetic  Resonance  ( NMR )  is  a  technique  which  probes  the  local 
environment  of  the  nuclei.  In  their  study  (15)  both  55Mn  and  27A1 
nuclei  were  examined  at  room  temperature.  One  of  the  conclusions  of 
this  study  supports  the  conclusion  of  the  Mossbauer  spectroscopy  (NGR) 
study  ( 13a )  . 

A  direct  observation  of  the  atomic  structure  of  the  Icosahedral  phase 
was  performed  by_Field_Ion_Microscopy  j.FlM_)_  (16).  The  -  nits  confirm 
the  long-range  Icosahedral  orientational  order  without  t-ms-  lational 
synsnetry.  The  basic  rotational  symmetries  were  observed  in  real 
space,  and  the  results  confirm  the  long-range  Icosahedral 
orientational  order  without  translational  synsnetry. 

In  addition,  a  large  number  of  apparent  discontinuities  were  found. 
Special  attention  was  paid  to  the  possibility  of  the  presence  of 
twins,  but  no  evidence  was  found  for  multiple  twinning,  and  no  twin 
boundaries  were  observed.  It  should  be  noted  that  the  FIM  studies 
were  performed  on  a  rapidly  solidified  ribbon  of  Al-12  at  %  Mn,  which 
is  by  nature  heavily  strained.  A  more  suitable  specimen  for  this 
study  would  be  a  rapidly  solidified  ribbon  of  the  Al-Mn-Si  ternary 
system  which  is  composed  of  elongated  strain  free  Icosahedral 
crystals . 

Based  on  the  FIM,  several  laboratories  have  now  developed  an  atom 
probe  which  is  potentially  an  important  tool  for  probing  the  chemical 
order  of  the  Icosahedral  phase,  and  for  evaluating  the  models  for  its 
atomic  structure. 

CONCLUSION 


Many  of  the  techniques  for  the  study  of  the  structure  of  solids  were 
applied  to  the  Icosahedral  phase.  The  long  range  order  was  determined 
by  X-ray,  neutron  and  electron  diffraction,  and  a  very  good  fit  was 
obtained  for  a  three  dimensional  quasiperiodic  arrangeawnt  of  atoms. 
The  local  order  was  probed  by  field  ion  microscopy,  nuclear  gamra-ray 
resonance,  nuclear  magnetic  resonance  and  extended  X-ray  absorption 
fine  structure.  At  this  stage  the  available  information  on  the  local 
order  should  provide  a  model  for  the  atomic  arrangement  both  of  the 
Icosahedral  phase  in  particular,  and  for  quasi-  periodic  crystals  in 
general.  Indeed,  such  models  are  proponed  in  these  proceedings,  and 
hopefully  an  agreed  upon  model  will  emerge  in  the  near  future. 
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COMMENTS  AFTER  THE  D.  SHECHTMAN  TALK  : 


M.  'J .  JARIC.  -  I  would  like  to  make  two  comments.  First,  X  would  like  ro 
emphasize  that  SHECHTMAN 's  model  of  randomly  packed  icosahedra 
shearing  edges,  faces,  or  vertices,  does,  in  fact,  exhibit  en  unusual 
kind  of  order  -  namely,  the  long-range  icosahedral  orientational 
order.  Second,  I  think  that  a  self consistent  definition  of  a  class  of 
(long-range  positionally)  ordered  structures  could  be  a  physical  one  : 
these  are  (infinite)  structures  which  give  rise  to  diffraction 
patterns  consisting  of  perfectly  sharp  bragg  peaks.  It  would  be,  then, 
an  experimental  task  to  recognize  this  ideal  in  real  structures. 
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materials  .  The  ribbons  were  then  gently  powdered  so  as  not  to  introduce 
excessive  strain  and  sieved  to  provide  the  neutron  diffraction  samples. 

III  -  EXPERIMENTAL  PROCEDURE 

Neutron  diffraction  measurements  were  performed  on  the  high  resolution  BT-1 
powder  diffractometer  /6/  at  the  N3S  20-MWatt  research  reactor.  The 
diffractometer  consists  of  5  separate  counters,  separated  by  20°,  that  can  move 
simultaneously  about  an  axis  at  the  sample  position.  The  5  counters  cover  an 
angular  range  from  0°  to  120°  and  for  our  experiments  covered  from  5°  to  120° 
scattering  angle  where  each  counter  moved  35°.  The  wave  length  used  for  our 
experiments  was  1 .5^6  A  provided  by  a  copper  monochromator.  An  oriented  graphite 
filter  was  used  in  the  beam  to  remove  higher  order  wave  lengths  contamination. 
The  resolution  of  the  diffractometer  varies  over  the  scattering  angle  from  0.31® 
at  5°  scattering  angle  to  a  minimum  of  0.22°  from  ^0°  to  70°  and  rises  to  O.fB® 
at  the  largest  scattering  angle. 

The  powder  to  be  studied  is  placed  in  a  thin  walled  vanadium  container  which  can 
be  fixed  at  the  sample  position.  Vanadium  was  used  since  it  does  not  contribute 
any  structure  to  the  diffraction  patterns.  Diffraction  data  was  taken  for  each 
powder  in  angular  steps  of  0.05°  and  for  a  given  monitor  count  of  the  incoming 
beam. 

IV  -  RESULTS 

Figure  1  shows  a  comparison  of  the  diffraction  data  versus  scattering  angle  for 
the  three  alloys  with  manganese  concentration  increasing  from  the  top  to  the 
bottom  graph.  The  lowest  concentration  manganese  alloy  when  melt  spun  under  the 
conditions  for  our  experiment  produces  the  icosahedral  phase  plus  the  f.c.c.  a- 
aluminum.  The  a-aluminura  phase  lines  in  our  pattern  are  narrower  than  the  peaks 
of  the  icosahedral  peaks  but  do  not  completely  reflect  the  instrumental 
resolution.  Their  width  could  arise  from  particle  size  broadening  or  possible 
hroadening  caused  by  non-uniform  concentration  of  Mn  ir.  the  a-aluminum  phase, 
ihe  icosahedral  lines  although  quite  sharp  compared  to  ether  samples  we  have 
examined  are  broadened  significantly  more  than  our  resolution  function.  One 
also  notes  an  oscillatory  background  under  the  peaks  of  the  icosahedral  phase. 
This  oscillatory  background  is  reminiscent  of  the  behavior  of  binary  alloy 
systems  showing  short  range  order  or  possibly  the  presence  of  an  amorphous  phase 
in  our  system.  No  amorphous  phase  was  detected  by  TZM.  An  examination  of  the 
background  behavior  at  the  smaller  angles  seems  to  indicate  a  pattern  similar  to 
short  range  order  since  the  curve  is  tending  upward  away  from  zero.  An 
experiment  is  being  designed  to  probe  these  smaller  angles  to  attempt  to  further 
clarify  this  point. 

One  sees  in  the  neutron  diffraction  a  slight  shift  in  peak  positions  versus  the 
manganese  concentration.  This  has  been  observed  by  both  electron  and  x-ray 
diffraction.  A/  With  increasing  manganese  content  there  is  a  change  in  the 
amount  of  other  phases  present;  a  aluminum  is  reduced  and  the  decagonal  phase 
appears.  (Neither  the  Al^Mn  nor  Al^Mn  phases  were  observed.)  These  changes  are 

observed  in  the  region  from  22°  to  i<50  in  two  thi>ta  but  can  be  seen  to  occur 
over  most  of  the  pattern.  We  are  attempting  to  prepare  a  sample  in  the  pure 
decagonal  phase  and  to  measure  its  neutron,  electron  and  x-ray  diffraction 
patterns  in  order  to  analyze  completely  the  diffraction  curves  with  mixed 
aperiodic  phases. 
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Fig.  1  -  The  uncorrected  neutron  scattering  intensity  of  the  3  alloys, 
Mn(27w/o),  AIMn  (30  w/o),  and  AIMn  (3*1  w/o),  from  top  to  bottom  respectively 
versus  the  scattering  angle,  28.  The  vertical  lines  are  the  predicted  positions 
of  the  Q  series  of  icosahedral  lines  for  all  even  values  of  N  between  2  and 


Figure  2  shows  a  comparison  of  a  powder  neutron  diffraction  pattern  of  the 
icosahedral  phase  with  the  X-ray  diffraction  pattern  /?/  for  a  similar  alloy. 
The  most  interesting  feature  is  the  difference  in  the  intensities  for  the  same 
icosahedral  peaks.  Furthermore,  many  peaks  are  observed  in  the  neutron 
diffraction  pattern  that  seem  to  be  absent  in  the  X-ray  pattern,  whereas,  some 
strong  X-ray  peaks  are  very  weak  in  the  neutron  pattern.  This  striking 
difference  in  intensity  comes  from  chemical  order  and  the  difference  in  phase  of 
the  neutron  scattering  amplitudes  of  the  manganese  atoms  which  have  a  negative 
neutron  scattering  length  relative  to  that  of  aluminum. 

In  the  preliminary  analysis  of  the  neutron  diffraction  data  for  the  275  alloy, 
we  obtained  the  integrated  intensity  under  all  icosahedral  peaks  positions 
calculated  from  the  index  method  we  discuss  later.  To  obtain  the  integrated 
intensity  we  did  not  attempt  at  this  point  to  fit  peak  profiles  witn  assumed 
mathematical  forms  but  subtracted  a  smooth  background  from  each  peak  and 
numerically  integrated  the  resultant  intensity.  For  some  positions  the 
intensity  differed  little  from  background.  These  were  listed  as  zero  intensity 
and  utilized  because  they  too  contain  information.  In  the  case  of  the  few 
overlapping  peaks,  where  the  separation  could  be  clearly  made,  we  adjusted  each 
peak's  intensity  to  give  it  a  smooth  symmetrical  curve  whose  summed  intensity 
agreed  with  that  observed.  Peaks  that  overlapped  with  those  of  other  phases 
were  not  measured.  In  order  to  compare  the  experimental  intensities  with  our 
model,  we  corrected  the  above  integrated  intensities  by  multiplying  it  by 

sin2e  cos 9 ,  the  applicable  Lorentz  factor  for  neutron  scattering,  where  0  is 
one-half  the  scattering  angle.  Finally  for  a  comparison  with  our  model  we 
divided  each  integrated  intensity  with  its  multiplicity  as  determined  by  the 
model  and  normalized  these  intensities  so  that  the  largest  one  would  be  unity. 
For  those  cases  where  several  nonequivalent  peaks  had  identical  0,  the  intensity 
was  divided  by  the  sum  of  the  multiplicities. 


Fig.  2  -  A  comparison  of  the  uncorrected  neutron  scattering  intensities,  top,  to 
the  scattering  intensities  for  x-rays  111  versus  Q  -  JJusinB/X  where  6  is  one- 
half  the  scattering  angle  and  X  is  the  wave  length  of  the  radiation  used.  The 
two  curves  in  the  neutron  data  near  Q*2.3  reflect  the  overlapping  range  of 
the  detectors 


V  -  DISCUSSION 


Two  aspects  of  the  neutron  diffraction  results  will  be  discussed;  peak  positions 
and  how  this  relates  to  the  icosahedral  symmetry,  and  intensity  and  hew  this 
confirms  the  emerging  model  of  a  structure  of  parallel  large  icosahedral  motif3 
stacked  aper iodically  along  their  3"fold  axes. 

VI  -  COORDINATE  SYSTEMS 

The  geometries  of  both  the  direct  and  reciprocal  spaces  in  three  dimensions  are 
conveniently  described  in  terms  of  a  Cartesian  coordinate  system  aligned  with 
one  set  of  three  mutually  orthogonal  2-fold  axes  which  occur  in  icosahedral 
symmetry. /8/  Each  of  the  three  components  of  a  quasilattice  or  reciprocal 
quasilattice  vector  will  be  described  in  terms  of  two  integers  e.g. ,  h  and  h’ , 
the  magnitude  of  the  component  is  h  +  h'x  where 

t  -  2Cos(H/5)  -  0  + "V5)/2  -  1  .618034. 


A  reflection  Q  is  indexed  with  six  integers  (h  +  h’x,  k  +  k’x,  1  +  l’x)  which 
will  be  written  (h/h*  k/k'  1/1').  Similarly  a  translation  vector  X  will  be 


(u/u*  u/u'  w/w').  These  six  integers  can  also  be  related  to  components 

6  ,  . 

.  x  )  of  a  six-dimensional  lattice  vector  in  Z  or  components  (n  n  . ..n,) 
6  •  £  u 
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of  a  reciprocal  lattice  vector  in  Z 


by  the  equations: 
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with  similar  equations  relating  u,u'...w'  to  the  x^. 

Equation  (1  )  can  be  considered  to  define  a  projection  or  Z  onto  a  three 

dimensional  plane.  It  will  be  convenient  to' let  the  (100000)  vectors  in  Z  and 

Z*6  have  a  unit  length.  Then  the  corresponding  vectors  in  R3  along  the  5-fold 
axis  (1/0  0/1  0/0)  will  have  a  projected  length  1  / .  We  therefore  normalize 
all  three-dimensional  vectors  by  a  factor  of  1 2(2  ♦  x).  In  our  previous 

paper  /8/  the  vector  (1/0  0/0  0/0)  in  R3  was  chosen  to  have  unit  length.  As  a 
result  of  the  redefinition  of  a  unit  length  formulas  for  lengths  in  this  paper 
differ  from  that  of  the  previous  paper  by  a  factor  of  2(2+x). 

Q  -  1  /S~2( 2  +  ~i)  ((h  +  h'x)i  +  (k  +  k’x)j  ♦  (1  +  l’x)k) 

X  -  1/12(2  +  x)  ((u  +  u'x)i  +  (k  +  k'x)j  +  (1  +  l’x)k) 

Note  that  X  (1/0  0/1  0/0)  -Z(  1  /0  0/1  0/0)  -  1/2.  These  vectors  which  were 

reciprocal  to  each  other  in  Z6  have  been  foreshortened  by  projection  and  are  no 
longer  reciprocal  to  each  other.  To  convert  these  to  dimensional  quantities  x 
and  q  we  will  use  the  six-dimensional  cubic  lattice  parameter  A  an  reciprocal 

* 

lattice  parameter  A  «1/A  and  let 


For  powder  diffraction  we  need  only  the  magnitudes  of  these  vectors.  Because 
t  ■  1  +  t  these  become 


Q  -  (N  +  Mt)/2(2  +  t)  X  -  (S  +  Tt)/2(2  +  t) 


and  where 


N  -  2Zn, 


S  -  2Zx, 


2  2  2  2  2 
+  h'  +  k  +  k'  +  1  +  1' 


2  2  2  2  2  2 
u  +  u'  +  V  +  v'  +  w  +  w' 


M  -  h'2+  k,2+  1'2+  2 ( hh ’  +  kk'  +  11')  T  -  u'2+  v'2+  w'2+  2(uu’+  vv*  +  ww' )  (5) 

Hence  N  and  S  are  always  even.  Furthermore  if  N(or  S)  is  divisible  by  M,  so  is 
M(or  T).  If  N(or  S)  is  of  the  form  Mm  +  2,  M(or  T)  is  of  the  form  Mm  +  1  .  In 
addition  we  have 


-N/t  <  M  <  Nt 


-  S/t  <  T  <  St. 


Equation  M  states  that  all  vectors  with  the  same  value  of  N  (or  S)  lie  on  the 
same  sphere  in  Z*6  (or  Z6).  The  distance  Qn  that  a  reflection  in  Z  6  is  from 
the  icosahedral  cut  plane  defined  by  equations  (1)  is  given  by 


CT  -  t(Nt  -  M)/2( 2  +  t) 
c 


Thus  the  largest  possible  value  of  M  for  a  given  N  will  have  the  smallest  Qc  and 
usually  the  greatest  intensity.  This  largest  value  of  M  is  called  Mq.  The 
corresponding  value  of  Q  is  called  and  is  a  one  parameter  sequence  of 


reflections. 


Gl  +  Mt|Nt/mJ 

2(2  +  t)  Q  =-  I  m  i  ( 

o  N  +  t(  1  +  ) 


N-Mn 


N-Mn+2 


where  j_xj  is  the  largest  integer  in  X  and  the  corresponding  values  of  Qc  are 


2(2+t)  Q‘ 


,T(r 


*6 

To  determine  multiplicities  of  the  reflections  all  nodes  in  Z  for  <  7  were 

grouped  by  their  length  Q  and  counted.  ^  ^ 

Table  1  lists  for  each  N,  the  largest  value  of  M,  the  indexes  in  Z  and  R  , 
the  multiplicities  and  the  values  of  Qq  and  This  same  table  can  be  used  to 

find  these  vectors  in  Z6  for  each  S  that  are  closest  to  being  parallel  to  the 
.  .  .3 


N  H  Coordinates  tn  R6  Coordinates  in  R3  Mult.  Q  ParaX.  q  perp. 
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Table  1 

The  main  distances  and/or  reflections  of  the  icosahedral  quasilattices  indexed 
according  to  N  and  M,  in  and  in  and  in  perpendicular  space. 


v:i  -  PEAK  POSITIONS 


In  figure  1  ,  we  compare  the  experimental  diffraction  curves  with  a  series  of 
vertical  lines  at  the  value  of  Qq  from  N  -  2  to  N  -  106  with  none  omitted  using 

A  -  0.65  nm.  Every  strong  peak  attributed  to  the  icosahedral  phase  coincides 
with  such  a  position.  A  few  weak  peaks  were  found  to  line  up  with  the  next 
lower  value  of  M,  M  *  M  -  Except  for  these  weak  peaks,  all  peaks  fall  at 

angles  determined  by  a  single  lattice  parameter  A=0.^6nm  and  one  other  number  N 

w'n.ch  takes  on  the  value  of  all  positive  even  integers.  In  a  sense  Q  (N)  is  the 

o 

icosahedral  equivalent  of  the  crystallographic  rules  which  determine  for  powder 
patterns  of  crystalline  materials  where  the  lines  are  to  be  found.  As  for 
crystalline  materials  it  has  a  beginning  at  N  -  2,  and  extends  monotonically  to 
higher  angles.  Strictly  any  list  of  all  peaks  must  involve  two  parameters  N  and 
M,  and  because  of  the  irrational  number  t,  Q  can  come  inf initesinally  close  to 
any  angle.  Only  the  main  Qq  sequence  has  high  intensity.  This  is  similar  to 

what  occurs  with  incommensurately  modulated  phases.  However,  the  strong  peaks 
from  an  incommensurately  modulated  phase  conform  to  a  periodic  lattice. 

VIII  -  INTENSITY 

We  find  that  the  intensity  of  neutron  diffraction  differs  significantly  from 
that  obtained  with  X-rays.  This  indicates  a  degree  of  chemical  ordering. 
The  intensities  are  not  consistent  witha  concept  of  a  simple  superlattice. 
since  there  is  no  systematic  intensity  difference  correlated  with  aby  parity 
rules  on  N.  In  any  ordering  scheme  the  N»1Ju+2  peaks  would  be  superlattice  peaks 
since  they  all  would  disappear  in  any  structure  that  belonged  to  an  icosahedral 
super  group.  As  in  ordinary  cryustallography ,  the  positions  of  the  reflections 
depend  entirely  on  the  six  dimensional  lattice  and  its  parameter,  only  the 
intensity  Is  affected  by  other  factors.  Because  the  icosahedral  phase  exists 
ever  a  range  of  compositions,  as  verified  by  the  shift  in  lattice  parameter, 
there  must  be  substitutions  of  manganese  and  aluminum  for  each  other  and  a 
certain  degree  of  chemical  disorder.  In  the  simplest  models,  identical  single 
scatterers  concentrated  on  a  quasilattice, /9, 10/  neutron  and  xr ay  diffraction 
would  give  similar  intensities,  which  would  be  a  known  function  only  of  Qc- 

This  kind  of  model  is  clearly  ruled  out.  More  complicated  models,  involving  a 
cuasilattice  decorated  with  identical  motifs  give  rise  to  a  factorable  structure 
factor  in  which  one  of  the  factors  gives  information  about  possible 
cuasilattices.  This  will  be  discussed  next. 

VIII  -  INFERENCE  ABOUT  THE  CUT  FUNCTION 

3 

A  quasiperiodic  arrangement  of  points  in  R  with  icosahedral  symmetry  is  easily 

obtained  by  projecting  all  lattice  nodes  in  within  a  symmetric  by  a 
neighborhood  of  the  icosahedral  plane  and  projecting  them  onto  the  icosahedral 
plane.  The  shape  of  the  neighborhood  must  conform  to  icosahedral  symmetry.  The 
size  and  shape  of  the  neighborhood  is  called  the  cut  function.  It  determines 
the  density  and  arrangement  of  quasilattice  nodes.  It  is  worth  noting  that  no 
points  in  the  quasilattice  will  be  separated  by  X  if  Xc  can  not  be  fit  into  the 

cut  function.  The  shape  and  size  of  the  cut  function  can  be  adjusted  to 
eliminate  specific  distances  in  the  quasilattice.  If  each  quasilattice  node  is 
decorated  with  the  same  atomic  motif  then  the  structure  factor  of  such  a 
structure  is  the  product  of  two  factors,  an  intrinsic  quasilattice  structure 
factor  and  one  derived  from  the  arrangement  of  atoms  in  the  motif.  The  former 
depends  entirely  on  the  cut  function.  The  separation  into  two  factors  is  an 
approximation.  In  the  ideal  structure  there  is  no  quasilattice  lattice  of 
points  which  have  strictly  identical  environments.  We  assume  an  approximate 
repetition  of  local  environment,  and  use  the  observed  intensity  to  make 
inferences  about  the  spacing  of  quasilattice  nodes  form  the  size  cf  the  cut 


Fig.  3  plots  intensity  vs.  t(Nt-M)  which  is  Qc  2(2+t).  The  quasilattice 
structure  factor  F  is  plotted  for  two  values  of  the  dimensionless  radius  of  a 
spherical  cut  function  R-1.1H  and  R-0.704.  The  former  has  a  node  in  Fig.  3  at 
1.65,  the  other  at  2.73*  The  larger  R  is  a  sphere  of  the  size  of  the  tri 
icoutahedroo  that  gives  the  quasilattice  of  the  3~dira  Penrose  tiling.  For  this 
size  the  principal  node  spacing  is  O.H6nm  along  the  5-fold  axis,  corresponding 
to  S-2,  T-1,  although  it  also  admits  a  few  of  the  short  diagonals  of  the  oblate 
tile,  S-6  T--3.  Because  of  the  many  intense  reflections  near  the  node  of  thi3 
structure  factor  we  conclude  that  this  is  not  the  size  of  the  cut  function,  that 
these  distances  are  not  quasilattice  vectors,  and  that  this  the  observed 
quasilattice  is  not  represented  by  the  3~diraensional  Penrose  filing.  The 
smaller  value  of  R  does  indeed  give  a  node  at  2.73  where  there  are  no  intense 
reflections.  This  size  cut  function  rules  out  S-2  and  S-4  and  all  S-6  except 
T-9,  which  becomes  the  shortest  distance  allowed  and  is  1.096nn  along  the  3-fold 
axis.  This  is  indeed  the  caliper  diameter  of  a  Mackay  icosahedron  /II/  along 
its  three-fold  axis  and  figures  in  the  model  proposed  by  Henley  and  Elser/1?./, 
Guyot  and  Audier,/13/  Ma,  Stern,  and  Eouldin  /1M/.  In  such  a  model  parallel 
Mackay  icosahedra  are  placed  on  nodes  separated  by  1.09Snn  along  their  common 
three-fold  axis  forming  distorted  octahedra  where  they  join.  No  nodes  are 
closer  than  this  distance.  In  another  paper  we  describe  the  construction  of  a 
cut-figure  which  will  give  a  quasilattice  with  these  properties. 
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